The traveling wave MRI in cylindrical Taylor- Couette flow: 
comparing wavelengths and speeds in theory and experiment 
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ABSTRACT 

We study experimentally the flow of a liquid metal confined between differentially rotating 
cylinders, in the presence of externally imposed axial and azimuthal magnetic fields. For in- 
creasingly large azimuthal fields a wave-like disturbance arises, traveling along the axis of the 
cylinders. The wavelengths and speeds of these structures, as well as the field strengths and 
rotation rates at which they arise, are broadly consistent with theoretical predictions of such a 
traveling wave magnetorotational instability. 

Subject headings: instabilities, magnetohydrodynamics, waves 



1. Introduction 

The magnetorotational instability (MRI) arises 
in a broad range of astrophysical problems, most 
importantly in accretion disks, where it is gener- 
ally accepted to be the source of the turbulence 
needed for material to spiral inward and accrete 
onto the central object (Balbus & Hawley 1991). 
Because of this crucial role that it plays in as- 
trophysics, there is considerable interest in trying 
to study the MRI in the laboratory (Rosner et 
al. 2004). One recent suggestion (Hollerbach & 
Rudiger 2005; Rudiger et al. 2005) involves ex- 
ternally imposing combined axial and azimuthal 



magnetic fields, which yields a new type of travel- 
ing wave MRI. In this letter we present experimen- 
tal evidence of these traveling waves, in accordance 
with the theory. 

The MRI is a mechanism whereby a differential 
rotation flow that satisfies the Rayleigh criterion, 
and is therefore hydrodynamically stable, may 
nevertheless be magnetohydrodynamically unsta- 
ble. The addition of a magnetic field allows angu- 
lar momentum to be transferred outward by the 
magnetic tension in the field lines, thereby bypass- 
ing the Rayleigh criterion, which relies on individ- 
ual fluid parcels conserving their angular momen- 
tum. The MRI is particularly relevant to Keple- 
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rian flows such as those found in accretion disks, 
where Q ~ r -3/2^ w hi c h WO uld indeed be stable 
according to the Rayleigh criterion. 

The simplest design for attempting to repro- 
duce the MRI in the lab is based on the familiar 
Taylor- Couette problem, consisting of the flow be- 
tween differentially rotating cylinders. While one 
cannot achieve precisely a Keplerian flow profile 
in this problem, by appropriately choosing the ro- 
tation rates of the inner and outer cylinders one 
can easily achieve a Ray leigh-st able profile, which 
is all that is really required. By suitably adjust- 
ing the strength of an externally imposed magnetic 
field one should then be able to destabilize the flow 
again, via the MRI. 

Unfortunately, the situation is not quite so sim- 
ple after all. If the imposed field is purely axial 
(Riidiger & Zhang 2001; Ji et al. 2001), the rel- 
evant parameter for the onset of the MRI turns 
out to be the magnetic Reynolds number Rm = 
Qirf /rj, which must exceed O(10). The hydrody- 
namic Reynolds number Re = Sl\r\jv then ex- 
ceeds O(10 6 ), due to the extremely small magnetic 
Prandtl numbers Pm = vjr\ of liquid metals [y is 
the viscosity, r] the magnetic diffusivity). Such 
large rotation rates can be achieved, but for in- 
creasingly large Reynolds numbers end-effects be- 
come increasingly important, and at Re > O(10 6 ) 
may well disrupt the experiment (Hollerbach & 
Fournier 2004). 

In contrast, if a combined axial and azimuthal 
field is imposed (Hollerbach & Riidiger 2005; 
Riidiger et al. 2005), the relevant parameter turns 
out to be i?e, which must only be as large as 
O(10 3 ) to obtain the MRI. These end-effects are 
therefore less severe. The solutions in this case are 
also somewhat different from those for purely axial 
imposed fields; one obtains much the same Tay- 
lor vortices as before, but the whole pattern now 
drifts along the length of the cylinders (see also 
Knobloch 1996). This unfortunately introduces 
significant end-effects of its own, familiar in other 
contexts, such as drifting dynamo waves (Tobias 
et al. 1998). Nevertheless, we will see that the 
experimental results in a bounded cylinder agree 
reasonably well with the theoretical results in an 
unbounded cylinder. 



2. Experimental Results 

The experimental apparatus consists of a cylin- 
drical annulus made of copper, with r\ = 4 cm, 
r Q = 8 cm, and height 40 cm. The top endplate 
is made of plexiglass, and is stationary; the bot- 
tom endplate is made of copper, and rotates with 
the outer cylinder. An axial magnetic field is im- 
posed by running a current, up to 200 A, through 
a series of coils surrounding the entire apparatus; 
an azimuthal magnetic field is imposed by run- 
ning a current, up to 8000 A, through a rod along 
the central axis. Field strengths of several hun- 
dred G can then be achieved, for both B z and 
B0. The fluid contained within the vessel is a 
GalnSn alloy, having density p = 6.4 g/cm 3 , vis- 
cosity v = 3.4 • 10 -3 cm 2 /s and magnetic diffu- 
sivity r] = 2.4 • 10 3 cm 2 /s, so Pm = 1.4 • 10 -6 . 
Measurements were made by two ultrasonic trans- 
ducers mounted on opposite sides of the top end- 
plate, 1.5 cm from the outer wall. These provided 
measurements of v z at these particular locations, 
along the entire 40 cm depth of the container. See 
also Stefani et al. (2006) for further details of the 
experimental setup. 

For the results presented here, the rotation 
rates of the inner and outer cylinders were fixed 
at fti = 0.377 s -1 and Q Q = 0.102 s~\ so 
Qo/fli = 0.27, within the Rayleigh-stable regime 
Q /Qi > (ri/r ) 2 = 0.25. The Reynolds num- 
ber Re = Vt\r\jv = 1775. The axial field B z was 
also fixed, at 77.2 G, corresponding to a Hartmann 
number Ha = B z r\j ppr\v = 12. The azimuthal 
field was then varied, between and 350 G (at 
n), corresponding to electric currents up to 7000 
A along the central rod. 

Figure 1 shows the temporal variation of the 
axial velocity v z , as a function of z along the con- 
tainer. The (depth-dependent) time-average has 
been subtracted, to remove the two Ekman vor- 
tices induced by the endplates (Kageyama et al. 
2004; Hollerbach & Fournier 2004). As B^ is in- 
creased, we see the gradual emergence of ever more 
coherent structures, drifting at speeds of 0.5-1 
mm/s. Focussing attention on z between 10 and 
20 cm, where the waves are most clearly defined, 
we obtain around 0.6 mm/s for 3 kA, and 0.8 
mm/s for 5 and 7 kA. The corresponding wave- 
lengths are 7-8 cm for 3 kA, and 5-6 cm for 5 and 
7 kA. We claim that these structures are precisely 
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Fig. 1. — Ultrasound measurements of the axial 
velocity v z , with the time-average removed. Time 
is plotted on the horizontal axis; height along the 
cylinder on the vertical axis. The current along 
the central rod varies from to 7 kA as indicated, 
corresponding to B^{r{) up to 350 G. 
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Fig. 2. — Contour plot of Re c as a function of 
B z and B^. On the horizontal axes B z is given 
both in G, and in nondimensional form, as Ha = 
0.1555 • B z j G. On the vertical axes B^(r\) is given 
both in G, and in terms of the required axial cur- 
rent I z . The numbers beside individual contour 
lines denote values from 1000 to 6000; the dashed 
contour line is 1775, the value used in the exper- 
iment. The five asterisks correspond to the five 
experimental runs in Fig. 1. 
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Fig. 3. — The left panel shows the wavelengths, 
in cm, the right panel shows the speeds, in mm/s, 
of the solutions at Re = 1775. As in Fig. 2, the 
dashed curves are the stability boundary at this 
particular Re. The three asterisks correspond to 
the three experimental runs at 3, 5 and 7 kA, 
where once again the wavelengths were approxi- 
mately 8, 6 and 6 cm, respectively, and the speeds 
0.6, 0.8 and 0.8 mm/s. 
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the expected traveling wave MRI. 

Note though that end-effects do indeed play an 
important role in Fig. 1, for example at the upper 
boundary, where the waves die away some 5-10 cm 
from the end. Another indication of the impor- 
tance of end-effects, and in particular the asym- 
metry between the two ends, can be seen in the 
A results: the asymmetry between top and bot- 
tom that is already visible even in this case must 
be caused entirely by the endplates, as a purely 
axial field does not distinguish between ±z. 

3. Theoretical Analysis 

Figure 2 shows the critical Reynolds number 
for the onset of this traveling wave MRI in an 
unbounded cylinder, as a function of the exter- 
nally imposed fields B z and B$\ these results were 
computed as in Hollerbach & Riidiger (2005) or 
Riidiger et al. (2005), with conducting boundary 
conditions. Provided B z > 30 G and > 150 
G, Reynolds numbers of O(10 3 ) are already suffi- 
cient. If B^ is less than 150 G, Re c gradually rises, 
until for B^ = we would have Re c > O(10 6 ), the 
familiar result from the analysis of purely axial 
fields. In contrast, if B z is less than 30 G, the in- 
stability simply ceases to exist; that is, Re c has an 
essentially vertical asymptote at this boundary. 

The five asterisks in Fig. 2 correspond to the 
five plots shown in Fig. 1, and the dotted con- 
tour line to the experimental value Re = 1775. 
We see therefore that the experimental runs with 
I z =0, 1 and 3 kA should be stable, whereas the 
ones with 5 and 7 kA should be unstable. This is 
broadly in agreement with Fig. 1, although there 
even the supposedly stable runs already show hints 
of traveling wave disturbances, particularly at 3 
kA. However, the waves are much more strongly 
developed for the 5 and 7 kA runs, as predicted 
by Fig. 2. 

Fig. 3 shows the wavelengths and speeds of the 
unstable solutions at Re = 1775. Wavelengths are 
in the range 10-20 cm, speeds around 1.1 mm/s 
(for B z = 77.2 G). The speeds agree rather well 
with Fig. 1; the experimental result that 5 and 7 
kA yield virtually the same speed, 0.8 mm/s, also 
nicely matches the theoretical prediction that the 
speed should be almost independent of B^ (pro- 
vided B§ is large enough to be in the unstable 
regime at all) . The wavelengths do not agree quite 



so well; the values in Fig. 1 are barely half those in 
Fig. 3. Presumably this is again due to end-effects; 
a 20 cm long wave certainly could not traverse a 40 
cm long container without experiencing significant 
end-effects. 

4. Conclusion 

In this letter we have presented experimental 
evidence for the existence of traveling wave dis- 
turbances in a liquid metal Taylor- Couette appa- 
ratus, and shown them to be in reasonable agree- 
ment with the theoretical predictions, particularly 
with regard to the wave speeds. Future experi- 
mental work will more thoroughly map out the en- 
tire B Z1 Be/,, Q{ and ^ Q parameter space, as well as 
explore the role that different axial boundary con- 
ditions, on both the flow (Kageyama et al. 2004) 
and the field (Liu et al. 2006), might play. Future 
numerical work will similarly consider the prob- 
lem in bounded cylinders, with different bound- 
ary conditions on the endplates. Only then will 
we fully understand what influence the endplates, 
and possible asymmetries between them, have on 
the MRI. 

This work was supported by the German Leib- 
niz Gemeinschaft, under program SAW. 
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